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We propose a new approach to the study of the inter-quark potential in Lattice Gauge Theories. 
Instead of looking at the expectation value of Polyakov loop correlators we study the modifi- 
cations induced in the chromoelectric flux by the presence of the Polyakov loops. In abelian 
LGTs, thanks to duality, this study can be performed in a very efficient way, allowing to reach 
high precision with a reasonable CPU cost. The major advantage of this numerical strategy is 
that it allows to eliminate the dominant effective string correction to the inter-quark potential (the 
Liischer term) thus giving an unique opportunity to test higher order corrections. Performing a set 
of simulations in the 3d gauge Ising model we were thus able to precisely identify and measure 
both the quartic and the sextic effective string corrections to the inter-quark potential. While the 
quartic term perfectly agrees with the Nambu-Goto one the sextic term is definitely different. Our 
result seems to disagree with the recent proof by Aharony and Karzbrun of the universality of the 
sextic correction. We discuss a few possible explanations of this disagreement. 
The numerical approach described above can also be applied to the study of Wilson loops. In this 
case, the numerical results are precise enough to test the two-loop prediction of the Nambu-Goto 
action. The two-loop NG result computed time ago by by Dietz and Filk is incompatible with the 
data; however, after correcting some mistakes in their expression, compatibility is restored. The 
viability of a first-order, operatorial description of the Wilson loop is also pointed out. 
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One of the most interesting recent results in the context of the effective string approach to 
LGT's is the proof of universality of the first few terms of the effective action. Universality up to 
the quartic order was proved a few years ago by Liischer and Weisz ||l|] and recently extended to 
the sixth order by Aharony and Karzbrun [^] (see also ^ for related results). These findings are 
based on a set of assumptions and it would be of great importance to test them with a numerical 
simulation. This is a very difficult task for at least two reasons: 

• in the standard "zero temperature" inter-quark potential, higher order corrections are propor- 
tional to higher powers in l/R{R being the inter-quark distance) and are thus visible only at 
very short distance, where the effective string picture breaks down and spurious effects (due 
for instance to boundary terms) and perturbative contributions become important; 

• the dominant string correction (the so called Liischer term) may shadow the sub-leading 
terms. 

In order to overcome these two problems we propose the following strategy. First, instead 
of working at zero temperature we shall study the inter-quark potential at finite temperature (just 
below the deconfinement transition). It is well known that in this regime the string corrections are 
proportional to R and act as a temperature dependent renormalization of the string tension. In this 
regime higher order corrections correspond to higher powers of T and can be observed much better 
than in the zero temperature limit. Second, in order to eliminate the dominant Liischer term we shall 
not measure directly the inter-quark potential, but the changes induced in the flux configuration by 
the presence of the Polyakov loops. We shall show below that as a consequence of this choice the 
Liischer term does not contribute. This makes this approach a very efficient tool to explore higher 
order corrections. 

The lattice operator which measures the flux through a plaquette p in presence of two Polyakov 
loops P, P' is: 

/pp't TJ \ 

{HP-,p,n) = ^-j^-{Up), (1) 

where Up is, as usual, the trace of the ordered product of Unk variables along the plaquette p. We 
shall be interested in the following to the mean flux density, defined as: 

where Np denotes the number of plaquettes of the lattice, R is the spatial distance between the 
Polyakov loops, namely, the inter-quark distance, and L the the number of lattice spacings in the 
time-like direction. In the following we will also introduce A^, to indicate the number of lattice 
spacings in the space-like directions, so that in (i = 2 + 1 we have Np = 3NjL. 

It is easy to see that if we define the partition function of the system in presence of the two 
Polyakov loops as 

Z(/?,L) = (P+(/?)P(0)) (3) 



then (<!>(/?, L)) can be written as: 
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where j8 is the coupling appearing in the Wilson action with respect to which the expectation values 
are being taken. If we neglect for the moment effective string corrections and keep only the area 
term in Z, i.e. Z{L,R) ~ e.^^^^ we find a linearly rising behaviour for (<!>(/?, L)): 

{^{R,L)) = aR (5) 

with an angular coefficient which in (2+1) dimensions is given by 

L do 1 do 

^ ^ N;,d^ ^ 3Njdp ^ ^ 

and does not depend on the finite temperature 1/L. 

The effective string corrections to Z{L,R) depend on the particular string action that we 
choose. They can be expanded in powers of the dimensionless quantity (oRL)^^: 

Z(L.;.)=e--^Z,.(l + ^ + ^ + ...), (7) 

where the indices in F4 and recall the fact that they are obtained from the quartic and sextic 
terms in the expansion of the effective string action respectively. 

The leading order of this expansion, namely Zi , corresponds to the partition function of a free 
boson in two dimensions. This term is universal and does not depend on the string tension o. Its 
dominant contribution in the large L limit is the well known "Liischer term" [Q]. The fact that Zi 
is scale invariant implies that it must also be j8 independent and thus disappears in eq. (^. As 
anticipated above, the effective string corrections to (<!>(/?, L)) start at the first sub-leading term, 
the "quartic" correction F4. In the Nambu-Goto case all the terms F„ are known [|l|, ||]. Inserting eq. 

in eq. @ we obtain the effective string corrections to (<!>(/?, L)). Since we shall be interested 
in the large R limit of these corrections, we find convenient to organize them via an expansion in 
powers of 1//?: 

{^{R,L)) = a(RAix)+B{x) + ^ + ..}j , (8) 

where x = In the Nambu-Goto case these functions can be evaluated explicitly ^ and turn 
out to be: 

(9) 

\/l—x 

While a depends on the details of the gauge theory and in particular it reflects the specific p 
dependence of the string tension o, the functions A{x),B{x),C{x), . . . encode the information on 
the effective string model. The particular form of these functions given in eq.s (^|llj) corresponds 
to the Nambu-Goto model; however, according to [|l|, ^], the first two orders in the perturbative 
expansion in power of 1/a (i.e. in powers of x) should be universal. In order to identify these 
terms, let us expand these functions in powers of x : 

A{x)=(l + ^ + ^ + ..) , (12) 
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^W = :r7 7 + ^ + - (13) 



oL \4 4 

CW = 3^(1 +.+ ...). (14) 

To test these corrections we performed a set of high precision simulations in the 3d gauge 
Ising model, using the same methods discussed in ^. We mapped via duality the Polyakov 
loops correlator into the partition function of a 3d Ising spin model in which we changed the sign 
of the coupling of all the links dual to the surface bordered by the two Polyakov loops. We then 
estimated (<!>(/?, L)) by simply evaluating the mean energy in presence of these frustrated links. 
Further details on the algorithm and on the results of the simulations can be found in Since 
duality plays a crucial role in this derivation, the approach discussed in this paper is particularly 
suited for abelian gauge theories; given enough computational power, however, it could be extended 
to non-abelian models. 

We chose to simulate the model at jS = 0.75180, for which both a and the deconfinement 
temperature are known with very high precision: a = 0.0105241(15) from [P] and 1 /T^ = = 8 



from [|10||. Using the above values for a and Tc and keeping into account the scaling correction 
as discussed in [11] we obtain for the prediction that a = 2.792 10^^ (see M\ for a detailed 



derivation). To test at the same time this prediction and the form of the effective string correction 
we selected two sets of values of L. The first set contains the values L = 16,20,24 and the variable 
X was tuned so that the sextic term x^/8 in eq. (12) is negligible with respect to the errors, but the 



quartic one, x^/8, is not. The second set instead comprises L= 10, 11, 12; in these cases also the 
sextic term is definitely not negligible with respect to the errors. 

For each value of L we fitted the data for <!>(/?, L) according to the law 

^{R,L) = a{L)R + b{L)+c{L)/R , (15) 

where the term c{L)/R was introduced only for the second set of values of L because it was always 
compatible with zero in the first set. In the following we shall concentrate on the values of a{L) 
which are the most precise and allow to perform a stringent test of the effective string prediction. 

The values of a{L) extracted from the fits are reported in tabjl] and plotted in fig. 1 . We analyzed 
these data in two steps. First we fitted the first three values of a{L) (those corresponding to L = 
16,20,24) with the law: 

a{L) = a{\ + Y^) . (16) 

We obtained a = 2.7918(17) 10"^ and y = 0.132(7) with a very good reduced x^- Both these 
values nicely agree with the predictions a = 2.792 10^^ and 7=1/8. 

We then fitted the whole set of data (i.e. including also L = 10, 11, 12) with the law: 

a{L) = a{l + yx^ + dx"") . (17) 

We found again a good reduced x^- Xr — 0-45 and the following best fit results: 

a = 2.796(5) 10"^ 7 = 0.127(25), 5 = -0.051(27) . 

The first two values agree again very well with the predictions but the coefficient of the sextic cor- 
rection, which should be 5 = 1/8, completely disagrees with them. This can be well appreciated 
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looking at fig.l, where we plotted the quartic correction (dashed curve), the sextic correction ac- 
cording to the Nambu-Goto effective action (dashed-dotted curve) and the curve corresponding to 
the best fit value of the parameter 8 (dotted curve). 



L 


a{L) 








10 


3.017(21) 


2.792 


3.137 


3.481 


11 


2.943(20) 


2.792 


3.027 


3.222 


12 


2.917(13) 


2.792 


2.959 


3.074 


16 


2.847(9) 


2.792 


2.845 


2.865 


20 


2.816(9) 


2.792 


2.814 


2.819 


24 


2.802(8) 


2.792 


2.802 


2.804 



Table 1: Values of the coefficient a{L) for various values of L. In the second column we list the results 
of the simulations extracted from the fits to eq.(15) In the following columns we report the prediction for 
a{L) = aA{x) at the zeroth, first and second order in the expansion in x 
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Figure 1: Plot of a A(L) as a function of jr = x- The continuous line is the prediction for a (i.e. without 
effective string corrections). The other curves correspond, from top to bottom to the whole Nambu-Goto 
prediction, eq.(9), the truncation at the sextic order, the truncation at the quartic order and finally the 
dashed dotted line corresponds to our best fit result (5 = —0.05 )for the sextic coefficient ( see the text). The 
points are the results of the simulations in the ltd gauge Ising model (see tab.l) 

This observation agrees with a set of similar results obtained in these last years in the 2x1 gauge 
Ising model by considering various physical observables, ranging from the inter-quark potential to 
the effective string width to the interface free energy of the dual spin Ising model [|l^, ^ [iT], |l^]. 
All these tests supported, even if only at a qualitative level, a value for the sextic correction different 
from the Nambu-Goto one and fully compatible with the one that we find here. 
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Our result is rather puzzling in view of the recent proof of the universality of the string ef- 
fective action up to this order in J = 3 A possible explanation could be that the observed 
deviation is only a lattice artifact due to unusually large scaling corrections and that the correct 
sextic contribution is recovered in the continuum limit. A careful analysis of the scaling behaviour 
of (<!>(/?, L)) is reported in ^ and seems to exclude this scenario: 5 shows a mild dependence on j8 
as the critical point is approached. It remains always negative and seems to slowly converge toward 
the value 5 = 0. 

Another possibility is that the deviation signals that 3d gauge Ising model does not admit a 
weakly coupled effective string description. This weak coupling limit is a basic assumption of 
^, ^ and amounts to ask that the partition function of the string describing a particular surface 
(say the cylindric surface connecting the two Polyakov loops) can be written as a sum of single 
string states propagating along the surface. The fact that this argument could be relevant for the 
3d gauge Ising model is also supported by the intuitive observation that a gauge theory based on 
the Z2 group is indeed the farthest possible choice with respect to the large limit of the SU(A^) 
gauge theory which is known to behave as a weakly coupled string theory. However if this is the 
reason behind the disagreement at the sixth order, it is not clear why the effective description works 
instead so well up to the quartic order. In order to better understand this issue we plan to perform 
the same analysis discussed here for other LGTs so as to gain some insight on the dependence of 
the 5 parameter on the gauge group. 

Another interesting application of our method is the study of the effective string corrections 



for Wilson loops. A long standing problem in this context is the fact that the Dietz and Filk [14] 



result for the quartic correction in the case of the Nambu-Goto action is manifestly incompatible 



with the Arvis spectrum [ ]15| ] for the open string. In particular, in the limit of very asymmetric 
Wilson loops, the Dietz and Filk result is one order of magnitude larger than (and opposite in 
sign to) the one suggested by the open string spectrum. To address this issue we performed a set 
of high precision simulations choosing the same values of j8 discussed above and found a value 
for the quartic term perfectly compatible with the Arvis spectrum while the Dietz and Filk result 
turned out to be excluded by more than ten standard deviations [|T6|]. Triggered by this result we 
went through the original Dietz and Filk calculation and identified the origin of the discrepancy. 
The correct result turns out to be compatible both with the Arvis spectrum and with the numerical 
estimates. For a rectangular Wilson loop of sizes R and T in a D + 2 dimensional LOT (i.e. with D 



transverse directions) it can be written in terms of Eisenstein functions as follows [16] 



F4 _ 1 /7r_Y 
oRT ~ oRT V24/ 



(18) 



The NG string corrections to the Wilson loop can also be computed along the lines used in 
and [ p^ in the case of Polyakov loops and interfaces, respectively. Replacing the Nambu-Goto 
second-order description with a first-order Polyakov one, the Wilson loop v.e.v. can be given an 
operatorial description: it corresponds to the propagation of a Dirichlet string whose end-points are 
attached to the spatial boundaries of the Wilson loop between two states which represent its emis- 



sion and re-absorption at times and T (such states were already considered in [18]). This way we 
can obtain a closed expression for the Wilson loop which can be expanded in powers of 1/ (aRT), 
which in the NG framework corresponds to the loop expansion; with respect to the NG treatment, it 
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is rather straightforward to extract higher loop contributions. Preliminary computations show that 



at two loop eq. ( |18| ) is reproduced up to the addition of a constant term in the square bracket; the 
precision of our numerical results does not allow to check if this constant term is supported by data 
or not. Such a (rather mild) discrepancy did not arise in the cases of Polyakov loops and interfaces, 



and it would certainly be very interesting to understand better its origin. In [ ]16[ ] we plan to develop 
the operatorial approach in a much more detailed way. 
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